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We calculate the cross section for diffractive hadron production in deep inelastic scattering off 
heavy nuclei in the framework of gluon saturation/color glass condensate. We analyze the kinematic 
region of the future Electron-Ion Colhder. We argue that coherent and incoherent diffractive channels 
are very sensitive to the structure of the nuclear matter at low x. This expresses itself in a character- 
istic dependence of the cross sections on rapidity and transverse momentum of the produced hadron 
and on the nuclear weight. We also discuss dependence on the scattering angle and argue that both 
coherent and incoherent cross sections may be within experimental reach at EIC. 

I, INTRODUCTION 

Diffraction is one of the most effective tools for investigating the structure of the nuclear matter at 
low values of Bjorken variable x. Its hallmark is large rapidity gaps (LRG) in rapidity distribution of the 
produced hadrons. At high energies, these gaps correspond to scattering processes mediated by exchange of 
a collective gluon state with vacuum quantum numbers, known as Pomeron. On the other hand, according 
to the Pomerantchuk theorem, high energy asymptotic of QCD is driven by the Pomeron exchange (see e.g. 
[1]). Hence, measurements of diffractive structure functions at HERA attracted a lot of interest. Indeed, 
diffractive physics at HERA yielded many exciting results that heralded the dawn of the new QCD regime 
of gluon saturation/color glass condensate (CGC) [2-11]. 

A possible launch of Electron Ion Collider (EIC) will open new avenues in studying the physics of 
diffraction in high energy nuclear physics. It will not only allow probing lower x and measure dependence of 
diffractive processes on nuclear weight, but also make possible studying less inclusive processes. One such 
process, diffractive hadron production in DIS is the subject of this paper. Our goal is to make predictions 
for DIS on a nucleus at the EIC kinematic region based on the CGC theory. We argue that diffractive hadron 
production is very sensitive to parameters of CGC and thus can be very effective instrument in extracting 
properties of the nuclear matter at low x. Gluon saturation effects on diffractive gluon production in DIS 
on proton at HERA have been discussed in [12-17, 23-25]. A concise discussion of the gluon saturation 
effects in semi-inclusive DIS on nuclei is given in [26-28]. 

This article is structured as follows. In Sec. II we review the formalism developed in our previous 
pubUcations [29-31], which allows to calculate coherent and incoherent diffractive gluon production in the 
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regime of coherent scattering » Ra, where Ic - l/(Mpx) is the coherence length in the nucleus rest 
frame. Coherent diffractive gluon production is the process y*+A^X + h + [LRG] + A. The corresponding 
cross section is given by Eqs. (l)-(3) and (6) below. For heavy nuclei A^^^ ~ l/aj » 1 and at high 
energies this type of diffractive process dominates over the incoherent diffraction, which is the process 
'y*+A^X + h + [LRG] +A with A* being excited nucleus. Nevertheless, at EIC energies, cross sections for 
coherent and incoherent diffraction processes are often comparable [31]. In pA collisions their dependences 
on gluon rapidity y and transverse momentum k and on atomic weight A are quite different. Therefore, as 
was pointed out in [31], it is important to separately measure the contributions of these diffractive processes. 
In Sec. m we calculate these contributions using the b-CGC model [35] for the color dipole scattering 
amphtude. As in [30] we characterize the nuclear effect using the nuclear modification factor (NMF) for 
diffractive processes defined in (1 1). The results of our numerical calculations are presented in Fig. 2. The 
most interesting features of the NMF's are (i) strong dependence of coherent diffractive NMF on gluon 
rapidity y (or xp); (ii) near independence of incoherent diffractive NMF on y and (iii) independence of both 
NMF's on the photon virtuality. This results are discussed in detail in Sec. III. 

Separation of coherent and incoherent diffractive contributions pose a great experimental challenge be- 
cause it requires measurements of very small scattering angles 6-2 sj-t/W^, where t is the moment trans- 
fer and W is the center-of-mass energy per nucleon of 7*A process. We address this problem in Sec. IV. 
Dependence of the coherent cross section on momentum transfer t is given by (24). It is seen that it de- 
creases as l/|?p at \t\ » ^/^A' where Ra is the nuclear radius. On the other hand, incoherent diffraction 
cross section decreases exponentially as e^'^'^p^^, but at much larger momentum transfers t > 1//?^ as seen 
in (36). The results of the calculation are plotted in Fig. 3. As expected coherent diffraction dominates at 
small momentum transfers -t while the incoherent one at large However, due to different functional 
form of ^-dependences, the two contributions become of the same order at about -t ~ Rp^ and remain 
comparable even at larger momentum transferes. The corresponding scattering angle for W - 100 GeV is 
w 0.13° and is very weakly dependent on the hadron transverse momentum, xp and photon virtuality Q^. 
It seems that such scattering angles are within the experimental reach and hopefully the two contribution 
can be separated. 
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II. DIFFRACTIVE GLUON PRODUCTION 



A. Dipole cross section 



Consider diffractive production of a gluon of transverse momentum k at rapidity y. Let the total rapidity 
interval be F = ln(l/x), where x - jW'^, is photon virtuality and W the center-of-mass energy of 'x*A^ 
scattering. Cross section for diffractive gluon production reads [18] 



da'^f^{r,x,k,y) 



cP-kdy 



/^*0'(Q,„) ^^^^^^ 



(1) 



where 



dcrl^ff{r,x,k,y) 



(2) 



d^kdy 

is the differential cross section for the diffractive gluon production by a qq dipole (a.k.a. onium) of trans- 
verse size r scattering off a nucleus. Eq. (1) generalizes the quasi-classical result derived in [19-21]. Other 
kinematic variables that are often used are yS and xp. They are defined as ln(l//3) = Y - y and ln(l/xp) = y, 
where Y - y is the rapidity interval between the photon and the produced gluon. We work in the approxi- 
mation ln( 1 /x) ~ 1, avln(l/y6) ~ 1. Diffractive production in the region yS < 1 was addressed in [12, 22]. 
We assume that the produced gluon is at the edge of the rapidity gap, so that the total rapidity gap in the 
process is y, see Fig. 1. 




FIG. 1 : One of the diagrams contributing to the diffractive production of a gluon with transverse momentum k and 
rapidity y. y is also the rapidity gap of the process. Unconnected f-channel gluons indicate all possible attachments to 
the i-channel gluons, quark and anti-quark. 
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Virtual photon light-cone wave-function appearing in (1) reads 



(D^ iQ,r,z) + ^l iQ,r,z) 

f 

2Ar,;^^42V(l-z)XV) 



(3) 




(4) 




(5) 



where = Q^z.{l - z) + m^p a{m = ^^Zy /(4:t), with Zf being electric charge of quark /. Subscripts L and T 
refer to the longitudinal and transverse polarizations respectively. 



We will consider two types of diffractive processes on nuclei - coherent and incoherent diffraction. 
Coherent diffraction is a process in which nucleus stays intact. This corresponds to elastic dipole scattering. 
At very high energies, such processes constitute half of the total dipole-nucleus cross section, another 
half being the inelastic processes. Therefore, contribution of coherent diffractive processes is expected to 
rise with energy. Unfortunately, experimental observation of coherent diffraction is challenging because it 
requires measurements at very small scattering angles, i.e. at very small momentum transfers \t\ ~ 1/^^. 
We discuss this in detail in Sec. IV. 

Another type of diffractive process is incoherent diffraction when the nucleus decays into colorless 
remnants. This process occurs at the nuclear edge where partial scattering amplitude at a given impact 
parameter is less than unity. Share of this contribution in the total inelastic cross section decreases with 
energy and with nuclear weight. Importance of incoherent diffraction stems from the fact that it measures 
fluctuations of the color glass condensate near its quasi-classical mean-field value. Typical momentum 
transfer in this case is \t\ ~ 1 /R^, i.e. determined by the inverse width of the diffuse region; it is much larger 
than in the case of coherent diffraction, which allows easier experimental study. In this section we discuss 
coherent and incoherent diffraction separately, assuming no experimental cuts on the minimal scattering 
angle. 

Cross section for coherent diffractive gluon production including the low-;c evolution was derived in 
[18, 29] and can be written as 



B. Coherent and incoherent diffraction 



do-cd(r,x,k,y) _ asCp 1 




d\' np{r, r',Y-y) |Icd(r', x, k, y, hf , 



(6) 
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where we introduced an auxiliary transverse vector 

X I - A^a(z - X, b, y) - NA{^ - y, b,y) + A^a(x - y, b, y) 

+Na(z - X, b, y) Na{z - y, b, y)} . (7) 

In [18, 29] we presented a detailed analytical and numerical analysis of the the coherent diffractive gluon 
production and discussed applications to pA scattering in [30]. Similarly, for incoherent diffraction [31] 

do-id(r, X, k, y) a^Cp 



(fikdy 2(2: 



-c- r r 

^Jd^bj d^r'nir,r',Y-y)pTAih)\Ii,,ir',x,k,y,bf 



(8) 



where 



/2 I z-x z-y\ 

xi^[l-NA{z-x,h,y)] [l-AfA(z-y,b,j)J [Np{z-x,0,y) + Np(z-y,0,y)] 

- [1 - NAix - y,h,y)] Npix - y, 0,y)} . (9) 

For numerical calculation we evaluate both vector functions Iia and led in the logarithmic approximation 
[18, 29-31]. Dipole density n(r, r',Y - y)d^r' in (8) is the number of daughter dipoles of size r' produced 
by a parent dipole of size r in the two-dimensional element of area d^r' at relative rapidity Y -y. In the 
diffusion approximation to the leading order BFKL equation [32, 33] it is given by: 



As discussed in detail in [30], nuclear modification factor Rab for coherent diff"ractive gluon production 
in the quasi-classical regime (i.e. without low-x evolution) is suppressed for large nuclei and large dipoles 
as Rqq+A ~ A^^^ exp{-r^(2^/4} (modulo logs) for dipole-nucleus scattering. Effect of quantum evolution is 
twofold. The larger is the rapidity of the produced gluon y, the stronger is the coherence effect that slows 
down growth of the diffractive qq +A cross section with energy as compared to the diffractive qq + p cross 
section. As a result, the nuclear modification factor gets an additional suppression in the y* fragmentation 
region (forward rapidity). On the other hand, at large 7 - y, the dipole density (10) in the virtual photon 
y* spreads to a wider range of sizes r' . Apparently, dipoles with sizes r' «; 21 Qs are not suppressed at 
all. This effect leads to relative enhancement of the nuclear modification factor in the backward versus 
forward rapidity. A quantitative study of diffractive hadron production requires numerical calculations that 
we discuss in the next section. 
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m. NUMERICAL CALCULATIONS 

A convenient way to express the nuclear effect on difitactive scattering is to introduce the nuclear mod- 
ification factor as a ratio of the difitactive cross sections in DIS on a nucleus per nucleon and on a proton 
[30]: 



cfikdy 

■'^cd/id - -r^ — ; • (11) 

^ iP-kdy 

Cross sections appearing in (11) are partonic cross sections (6) and (8) convoluted with the LO pion frag- 
mentation function given in [34] . 

We performed numerical calculations with b-CGC model of the dipole scattering amphtude A^^ [35], 
albeit with a modification: we treat nuclear and proton profiles as step-functions; the saturation scales are 
assumed to scale with A Q^^ A^l^ . The advantage of this model is that (i) its form compUes with the 
known analytical approximations to the BK equation and (ii) its parameters are fitted to the low x DIS data. 
The expUcit form of the scattering amplitude is given by 

N{Y,0,y)= \ ^ I (12) 

[ 1 - exp[-a \n\brQs)\ , rQ, > 2 , 

where is the the quark saturation scale related to the gluon saturation scale 2^ - which we have referred 
to simply as the 'saturation scale' throughout the paper - by = (4/9) 2^- Its functional form is 

Ql^A^l^x'^e^^GQN^, (13) 

The anomalous dimension is 

r = y, + ^ln(-^]. (14) 

Parameters js - 0.628 and k - 9.9 follow from the BFKL dynamics [36], while M) - 0.7, xo = 3 ■ 10"^ 
and A - 0.28 are fitted to the DIS data. Constants a and b are uniquely fixed from by the requirement of 
continuity of the amplitude and its first derivative. 

Our results are presented in Fig. 2 which exhibits dependence of the nuclear modification factor for 
coherent (left column) and incoherent (right column) hadron production on transverse momentum k. We 
assumed that the center-of-mass energy of the 7* A collision isW - 100 GeV per proton, which corresponds 
to the total rapidity interval Y - 9.2. 

In Fig. 2 (a,b) we show variation of the nuclear modification factor with the nuclear weight. We observe 
that ^cd increases with A. This is a signature behavior of higher twist eff'ects and, in particular, coherent 
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W=100 GeV, Q^=l Ge\\ xp=0M67 



W=100 GeV, e^=l GeV^, xp=0.0067 




(a) 

W=100 GeV, Q^=l GeV^ A=197 
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(c) 

W=100 GeV, ji;p=0.0067, A=197 
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W=100 GeV, Q^=l GeV^ A=197 
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W=100 GeV, jf =0.0067, A=197 
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FIG. 2: Nuclear modification factors for coherent (left column) and incoherent (right column) difFractive hadron 
production at W = 100 GeV as a function of the hadron transverse momentum k±. Shown are dependences on: (a),(b) 
atomic number A, (c),(d) hadron rapidity y and (e),(f) photon virtuality Q^. 



diffraction. In view of the discussion at the end of the previous section, we infer that the effective dipole size 
r' produced in the dipole evolution is r' «; 2/Qs, for otherwise the cross section would decrease for heavier 
nuclei. As one can see in Fig. 2 (e,f), NMF has no significant dependence, and hence no r dependence 
as well. Therefore, even at higher y, where evolution effects in the nucleus as well as lack of evolution in y* 
could have produced suppression of Red with A, no such suppression is observed. We checked this statement 
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up to the most forward direction allowed by our model j3 - OA. decreases with A already at midrapidity 
y - 5 because the general property of incoherent diffraction is that it vanishes in the Umit A — > oo when all 
partial amplitudes turn black. 

Rapidity dependence is displayed in Fig. 2 (c,d)- Red rapidly decreases in the forward direction, which is 
a cumulative effect of evolution in the nucleus and in the virtual photon, whereas Rid is essentially rapidity 
independent. This effect has already been noticed by one of us in pA case [31]. It arises because of different 
physical origins of the two diffractive processes. Coherent diffraction corresponds to elastic scattering of 
a color dipole on a nucleus, whereas incoherent diffraction is a part of inelastic scattering that originates 
from the nuclear periphery due to variation of the nuclear density with impact parameter. At low x central 
impact parameters of a heavy nucleus are black for a typical dipole. Therefore, scattering amplitude of 
dipole on a heavy nucleus is very different from an incoherent superstition of dipole-nucleon scattering 
amplitudes, hence strong variation of the nuclear modification factor with energy/rapidity. On the other 
hand, incoherent diffraction is non-zero only in the range of impact parameters comparable with the proton 
radius. Therefore, energy/rapidity dependence of dipole-nucleus and dipole-proton cross section is similar, 
though the geometry is quite different. 

Finally, Fig. 2 (e,f) exhibits dependence on photon virtuahty Q^, or perhaps better to say no dependence 
at all. This can be interpreted as insensitivity of the diffractive cross sections to the size of the parent dipole 
r. Indeed, as explained in [29], at k± » Qs, Q diffractive spectra depend only on k_L. For example, cross 
section for coherent diffractive gluon production in the asymptotic kinematic region Qs l/r <^ k reads 
(in the double-logarithmic approximation) 

^ = ^l££lAN\l/k,b,y) ^ 2V2a.(r-y)in(.fc)_ (15) 

d^kdy n^l^k^ ^ {Id^iY -y)\n{rk)yi^ 

Clearly, r-dependence cancels out of the nuclear modification factor. Notice, however, that the EIC kine- 
matic region can hardly be classified as asymptotic, and one would expect large corrections to (15). In 
fact, it is known that corrections to the double-logarithmic approximation are phenomenologically signifi- 
cant (see e.g. [16, 37]). However, our numerical calculations imply that they cancel in this particular case. 
Unfortunately, we are not able to extend this analysis to higher Q^'s without transgressing the region of 
applicability of our model. It would be interesting to analytically investigate the origin of this cancelation. 



IV. T-DEPENDENCE 



In this section we consider dependence of different diffraction channels on momentum transfer t. t- 
dependence translates into dependence on the scattering angle 6. While the dominant contribution to the 
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difitactive cross sections stems from scattering at small angles, only angles larger than some cutoff angle 
are experimentally accessible. In this section, we would Uke to investigate whether separation of coherent 
and incoherent contributions is experimentally feasible at ETC. 



A. Coherent diffraction 



Consider dipole-nucleus elastic scattering amplitude T'<'^'^^(s, b, {ha}), where b is the dipole impact pa- 
rameter and bfl's are positions of nucleons in the nucleus. Average over the nucleon positions will be 
denoted as (t'^^{s, b)^. Cross section for elastic dipole scattering is 



qq+A 



= J d^b \{T'^^^\s,h))f . (16) 



In this representation, (16) is also the coherent diffraction cross section. Fourier image of the dipole-nucleus 
elastic scattering amplitude carries information about the transferred momentum A(t - -A^): 

A)) = 2 J d% b)) e^ ^ . (17) 

If only two-body forces are taken into account in the scattering amplitude, which amounts to neglecting cor- 
relations between nucleons, then we can express the scattering amplitude on a nucleus through the scattering 
amphtudes on individual nucleons as 

A 

ri^^\s, b, {ba}) = 1 - ]~[ (i - r^^+^(5, b - bj) . (18) 

a=\ 

In this approximation, averaging can be performed as 

{...)^W dX dzpAiha, z)... = Y\ dXpTAiK) ... (19) 

„ 1 *J *J—oo „ 1 %J 

a=l a=l 

where pA(b, z) is the nuclear density at a given point in the nucleus and p is its average over the nucleus 
volume. 

Impact parameter profile of the dipole-nucleon amplitude is traditionally parameterized as 

n^+^(., b) = ^-^Z'^'is) -^e-'^l^l , (20) 
2 nRj 

where we neglected a small imaginary part of T'^^'^^is, b). In a heavy nucleus of radius Ra » Rp, nucleon 
can be approximated by the delta function in impact parameter space. Thus, 

J d^V^^^is, b - ha)pTAiha) ^ pTi^^^is, 0) pTA(h) . (21) 
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Using (19),(20),(21) in (18) we derive for heavy nuclei 

, b)) = 1 - 5'-tr"(^)p7'A w (22) 



Finally, substituting (22) into (17) and (16) we find 



dt 



._i.|2j j2^(l-e4-r"(W)),'b.. 



(23) 



To estimate the ^-dependence of the coherent cross section we can use a simple model for the b- 
distribution. Denote \(r^^^ {s)pTA{b) - O.S{b) and let the profile function S{b) be given by the step 
function 5 (Z?) = 6{Ra - b). Neglecting contribution of the diffuse region at the nucleus edge is a reasonable 

approximation in the case of coherent diffraction because the main contribution stems from b < Ra impact 
parameters. Substituting into (23) and (16) we get the well-known result (see e.g. [1]) 



dt a-T^ \t\ 

ca 



(24) 



Because (24) does not depend on Q this formula also gives ^-dependence of the diffractive coherent gluon 
production: 

do-ll^{Q\x,k,y) _ jI{Ra ^Tt) da-l'^\Q\x,k,y) 
d^kdydt ~ \t\ d^kdy 

B. Incoherent diffraction 

Coherent diffraction includes only events in which nucleus stays intact. However, generally the nucleus 
can be excited and subsequently decays into colorless remnants. Total diffractive cross section of this 
process is given by 



= J d^b i^Y^~^^\s,hfY (26) 



The difference between (26) and (16) measures dispersion of the scattering amphtude in the impact param- 
eter space. The corresponding physical process is a part of inelastic cross section and is called incoherent 
diffraction: 

trff = J d^b i^Y'i'^*\sMt) - |(r^^+^(^,b))|^ . (27) 

Clearly, the incoherent diffraction stems from the region near the nucleus edge ('diffuse region'). Indeed, 
dXb <^ Ra all partial dipole-nucleon amplitudes are close to the black disk limit, while at ^7 » /?a they all 
vanish. 
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To derive the ^-dependence of the incoherent diffraction cross section we define similarly to (17) 



(28) 



Then (26) reads 
do-dif 



dt 



(29) 
it. 



+g-i:<,<r««*^(i,b-b<,)>+i:„(r«+JV(^^b/_bj)_^r«5+'V(j,b-b„)r««+JV(i,b'-ba))j _ 



Upon subtracting the coherent difiraction part 
dcTcd 1 



dt An 
we end up with 

JcTid 1 



J d'b j jWMi'-b')(i_g-i:.<r«-^(..b-b„))j^i_^-i:, 



(30) 



(31) 



^ LJd2^J ^2^,^;A.(b-b')[i_^-i:.<r--A'(.,b-b„)n?-^(.,b'-bj)J 
xe-2a[<r«+"(i,b-bj)+<rw+"(i,b'-b„))-<r«9+'V(^,b-b„)r«+"(i,b'-b„))] _ 



(32) 



Since elastic qqN cross section is small compared with the inelastic one (as it contains extra [31]), we 



expand (32): 



dt 



do-id ^ ]_ J^if^ J ^2^/giA (b-b')g-Za[<r«^''(^,b-b„))+<r««+«(5,b'-b„))] 
X J] {n^+^(^, b - b„)n^+^(^, b' - ha)) 

a 

= 1 J j2^7g'^^-''^-*(^)r*'^''(^'0)n^^^(b-b« 



Because |b - ba| ~ Rp ^ ba ~ Ra, (34) becomes 



JcTid 1 



-2prA(6„)r??-'V(i,0) 



dt An , 

Finally, using (20) we derive the desired result 



1/ ^^^''^ 



bp^'^+A'^jj-j 



Jo-id 1 o-tor^('^) _-i 



prA(ba) 



prA(ba) . 



e-'^'^l J d% e-'^^^^''-^^''''^^%TAha) = ^.-^"1^? . (36) 



(33) 
(34) 

(35) 



dt An 2 

As in the case of coherent difiraction, ^-dependence of the cross section for incoherent diffraction is inde- 
pendent of other kinematic variables and therefore (36) describes also ^-dependence of incoherent diffractive 
hadron production: 



dcrj/(Q^,x,k,y) _ dcrJ^^(Q'^,x,k,y) Rj ^_i^,^2 



d^k dy dt 



d^kdy 



(37) 
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61° 




-t(GeV^) 

(a) (b) 

FIG. 3: (a) f-dependence of difFractive coherent (solid line) and incoherent (broken line) hadron production, (b) 
Contour plot of the scattering angle 0' at which coherent and incoherent difFractive hadron production cross sections 
are equal at given k and y = ln(\/xp). 

f-dependence of diffractive coherent and incoherent hadron production is shown in Fig. 3. In Fig. 3(a) 
we compare the coherent and incoherent difFractive production as a function of the transferred momentum t 
for a particular choice of the collision kinematics. We observe the following general features: (i) At small |f| 
coherent diffraction dominates over incoherent one by several orders of magnitude; (ii) At \t\ ~ coherent 
and incoherent contributions coincide, we will call the corresponding scattering angle 9'; a.t\t\ > R^^ the 
two contributions are similar on average; (in) Features (i) and (ii) hold for a wide range of parameters. In 
particular, 9' is nearly constant 9' w 0.13" as is seen in Fig. 3(b). Scattering angles of this size may be 
experimentally accessible at EIC. 



V. SUMMARY 



In this paper we discussed coherent and incoherent diffractive gluon production in DIS off heavy nuclei 
in the proposed kinematic region of Electron Ion Collider. Our approach is based on the dipole model 
introduced in [38]. It allows representing cross sections for high energy hadronic scattering as a convolution 
of hadronic light-cone wave-functions with the multipole scattering amplitudes. In our case, virtual photon 
wave function is determined by the perturbative QED and is given by (3). Dipole-nucleus interaction can 
in turn be represented as a product of dipole density (10) in transverse coordinate space, satisfying the 
BFKL equation [32, 33], and the dipole-nucleus forward elastic scattering amplitude as displayed in (6),(8), 
satisfying QCD evolution equations in the low x region [6, 39]. These formulas are derived in the leading 
logarithmic approximation asln{l/x) ~ 1, asln{\/j3) ~ 1, which defines the kinematic region where the 
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results of our calculations are applicable. Note, that hard perturbative factorization is generally broken 
at low X, because scattering in this region is characterized by small longitudinal momentum transfer (see 
e.g. [40]). At moderate x and large (2^, our formulas reduce to the leading order hard perturbative QCD 
expressions that can be cast in the factorized form using the diffractive parton distributions [41^4]. 

The main results of our calculations are displayed in Fig. 2 and Fig. 3. We found that nuclear modifica- 
tion factor strongly varies with nuclear weight, and the functional dependence on A is qualitatively different 
for coherent and incoherent processes. Similarly to diffractive hadron production in pA collisions [31], 
nuclear effects in coherent diffractive DIS is strongly dependent on rapidity of produced hadron, whereas 
they are almost absent in the case of incoherent diffraction. We also made a pecuUar observation that the 
nuclear modification factor for both diffractive channels is essentially independent of the photon virtuaUty 
in the region \<Q^ <25 GeV^. Finally, our study of non-forward diffractive hadron production indicates 
feasibihty of experimentally separation of coherent and incoherent diffractive contributions at EIC. 
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